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Abstract. Electrodynamics for self-interacting scalar fields in spatially flat Friedmann- 

00 ; 

Robertson- Walker space-times is studied. The corresponding one-loop field equation for 

: 

the expectation value of the complex scalar field in the conformal vacuum is derived. For 
5h . exponentially expanding universes, the equations for the Bogoliubov coefficients describing 

bX). 

.^J ' the coupling of the scalar field to gravity are solved numerically. They yield a non-local 
correction to the Coleman- Weinberg effective potential which does not modify the pat- 
tern of minima found in static de Sitter space. Such a correction contains a dissipative 
term which, accounting for the decay of the classical configuration in scalar field quanta, 
may be relevant for the reheating stage. The physical meaning of the non-local term in 
the semiclassical field equation is investigated by evaluating this contribution for various 
background field configurations. 
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1. Introduction 

In a recent series of papers [1-3], some of the authors have studied the one- loop effective 
potential for grand unified theories in de Sitter space. Our main results were a better 
understanding of the symmetry-breaking pattern first found in [4], a numerical approach 
to small perturbations of de Sitter cosmologies [2], and the analysis of SO(10) GUT theories 
in de Sitter cosmologies [3]. However, a constant Higgs field, with de Sitter four-space as 
a background in the corresponding one-loop effective potential, is only a mathematical 
idealization. A more realistic description of the early universe is instead obtained on 
considering a dynamical space-time such as the one occurring in Friedmann-Robertson- 
Walker (hereafter referred to as FRW) models. Indeed, our early work appearing in [2] 
tried to study the case of varying Higgs field by introducing a function of the Euclidean-time 
coordinate which reduces to the 4-sphere radius of de Sitter in the limit of constant Higgs 
field. Although the approximations made in [2] were legitimate for numerical purposes, 
the gravitational part of the action, and the one-loop effective potential, were not actually 
appropriate for studying a dynamical cosmological model. 

Thus, relying on the work in [5-6], this paper studies the first step towards the com- 
pletion of our programme, i.e. scalar electrodynamics with a self-interaction term for the 
complex scalar field. The semiclassical field equations in time-variable backgrounds con- 
tain non-local terms, which describe the coupling of the scalar field to the gravitational 
background. The analysis of these equations is relevant for the reheating mechanism in 
inflationary cosmology and for the dynamics of dissipation via particle production [7] . Our 
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analysis deals with Lorentzian space-time manifolds with FRW symmetries, and does not 
rely on zeta-function regularization. Hence our geometric framework is substantially dif- 
ferent from the Riemannian four-manifolds studied in [1-4] , where the metric was positive- 
definite. Section 2 derives the field equations for a class of cosmological models where 
scalar electrodynamics is studied in FRW universes. The Coleman- Weinberg potential, 
and its correction deriving from the non-local term involving the Bogoliubov coefficients 
for the coupling of the scalar field to gravity, are obtained numerically in section 3. The 
semiclassical field equation and the physical meaning of such non-local term are studied 
in detail in section 4. Concluding remarks are presented in section 5, and relevant details 
are given in the appendix. 

2. Model and field equations 

For the reasons described in the introduction, we consider a complex scalar field (p with 
a mass and a self-interaction term, coupled to the electromagnetic potential in curved 
space-time. Hence the action functional is 




boundary terms 



(2.1) 



where (cf [6]) 



C 



(2.2) 
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With a standard notation, V is the Levi-Civita connection on the background space- 
time, £ is a dimensionless parameter, R denotes the trace of the Ricci tensor, and F^ v = 
V V A^ — V^Ajy is the electromagnetic-field tensor. Boundary terms are necessary to obtain 
a well-posed variational problem, and their form is obtained after integration by parts 
in the volume integral in (2.1). Covariant derivatives of the scalar field are here used to 
achieve a uniform notation (cf [6] and [8]). 

Following [5], we now split the complex scalar field <fi as the sum of a variable, real- 
valued background field C , and of a complex-valued fluctuation <p, i.e. 



The conformal vacuum [9] is here chosen, and the quantum fluctuation <p has vanishing 
expectation value in such a state, < <p >= 0, so that < >= C . The field equations for 
Ap, 4> c and <p are obtained by setting to zero the corresponding functional derivatives of the 
action. As far as the gauge potential A^ is concerned, it is convenient to impose the Lorentz 
gauge VM^ = 0. At this stage, to quantize the theory, one can follow the Gupta-Bleuler 
method, or the Faddeev-Popov procedure, or to eliminate the residual gauge freedom by 
imposing the relativistic gauge condition proposed by Ford in [6], i.e. X^A^ = 0, where 
X is a timelike vector field. Such a field admits a natural form in FRW space-times (see 
below). Thus, on defining the operator [~] = g^ u V ^V v , the resulting form of the field 
equations is (cf [6]) 



(j) = (j) c + tp. 



(2.3) 




(2.4) 
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□ + m 2 + £R + ^-Xcfl - e 2 < A^A* > + -*-A < cp 2 + (v? f ) 2 + 4 W f > <f> c 



ie < A^V^ip - ip ] ) > 



(2.5) 



□ + m 2 + £i? + |A^1 ^ + J_ A 2 = o. 



(2.6) 



Note that equations (2.4)-(2.6) have been obtained by retaining in the action (2.1) only 
terms quadratic in the fluctuations <p and A^, and setting to zero all terms involving 
V v A v and its covariant derivatives in the field equations. The latter condition is sufficient 
to derive (2.4). Moreover, we require that V M C should be proportional to AT M [6]. Note 
also that the contribution of 



has been neglected in the course of deriving (2.6), since equation (2.5) implies that such 
a contribution is of second order in the quantum fluctuations. By taking the complex 
conjugate of equation (2.6), and defining <p = (<pi + i<p2)/V2, Xi = A/2 and A2 = A/6, the 
addition and subtraction of the resulting equations leads to decoupled equations for the 
real and imaginary parts of <p, i.e. 



Moreover, by virtue of our particular gauge conditions, equation (2.5) takes the form (see 
the appendix) 



□ +m 2 +tR+— X& C 




(2.7) 
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Interestingly, the effects of quantum fluctuations in equation (2.8) reduce to a linear super- 
position of the self-interaction term studied in [5] and of the electromagnetic term studied 
in [6], without any coupling. The term < A^A^ > is evaluated on considering the integral 
equation equivalent to (2.4), as shown in [6] and in the appendix. Such an analysis proves 
that, in the case of a spatially flat FRW background, the renormalized expectation value 
of A^A^ in the conformal vacuum can be written as [6] 

< A^ >= B X^X v R tw (2.9) 

where B = 9.682 ■ 10 -3 . In particular, in a de Sitter universe, < A^A^ >= 12BH 2 , where 
H is the Hubble parameter [6] . 

Following [5], we can now write the set of equations leading to the solution of (2.7) 
and (2.8) in spatially flat FRW backgrounds. Without making any approximation, if one 
defines (a being the cosmic scale factor) 



(2-10) 



and, for all j = 1,2 



fl5 ifc = fc 2 + o a (r) 



m 2 +^-l)R(r) + l\ j( p 2 c (T) 



(2.11) 



the one-loop field equations resulting from (2.7) and (2.8) are (cf (3.32) and (3.33) in [5]) 

l_rf^ £ + ^_rfa# £ _ e2< A AM>0c+ 3Veff 
a 2 dr 2 a 3 dr dr c d(j) c 

(b 2 1 f°° r i 

+ E 2 A ^ I dk fc2 °^ + Re z >A = (2 - 12) 
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d 



>j,k = I -j- log %fc Re z jik 




(2.13) 




( 



— log % fc J (s i)fc + - J - 2i Q jik ^ >fc . 



(2.14) 



The form of dV e s/d(f) c is given in the appendix. Equations (2.13) and (2.14) are necessary 
to find solutions of (2.7) by using Fourier-transform techniques and a suitable change of 
coordinates, as shown in [5]. The initial conditions for Sj^ and Zj jk are the ones appropriate 
for the choice of conformal vacuum, i.e. [5] 



Of course, the values of m, £ and Xj should be now regarded as the renormalized values of 
such parameters [5]. 

Note that the second line of equation (2.12) is a typical non-local term, resulting from 
the self-interaction of the scalar field (see (2.7)) and from its coupling to the geometric 
background. Moreover, a dissipative term exists which is part of the non-local correction 
to the Coleman- Weinberg potential, and it is due to the decay processes of <j) c in scalar field 
quanta. By means of such decays, energy is transferred from 4> c to the relativistic degrees 
of freedom ip (i.e. radiation). As it is well known, if this release of energy is sufficiently 
strong, radiation becomes dominant and hence the inflationary phase ends. This leads to 
the reheating stage, which is as important as the exponential expansion for the dynamics 
of the early universe. 



8 j>k (T = 0) = 



z i>fc ( T = 0) = 0. 



(2.15) 
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3. Numerical evaluation of the one-loop effective potential 



In this section, for the physically relevant case of an exponentially expanding FRW uni- 
verse, we compute the non-local term on the second line of (2.12). After integration of such 
a term with respect to C , the resulting expression is compared with the Coleman- Weinberg 
potential V e Q. 

The numerical analysis has been performed for a conformally invariant scalar field 
(hence £ = 1/6 and m = 0), by solving equations (2.13) and (2.14) with the help of the NAG 
routine D02BAF, with initial conditions (2.15). We have fixed A = 10~ 2 , which ensures the 
reliability of the perturbative approach, and the Hubble parameter H = 10 _1 Mpl, which 
is an intermediate choice between a chaotic model (H ~ M PL ) and a GUT inflationary 
phase (H ~ 10" 4 M PL ). 

In the case of exponential expansion in the time variable t, which implies (see (2.10)) 



and for a classical field configuration independent of r, C = c o, we express the second 
line of (2.12) as a function of r and 4> c q. In figure 1 we plot the integral with respect to 
4> c q of the above quantity when the conformal time r varies between and 1/H, which 
corresponds to a large e-fold number, cO G [0, 10 M PL ], and a(0) has been set to 1. This 
configuration for o(r) and fixed c o corresponds to an exact de Sitter space, but unlike [1- 
4] , with a Lorentzian signature. The non-local term describes an energy exchange between 
the gravitational field and the field <f> (see (2.3)). Such an exchange may lead to dissipative 
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or non-dissipative processes in the early universe, depending on the sign of the non-local 
term, as it will be clear from the analysis in section 4. For a constant background field <p c , 
only a scale factor varying in time yields the above non-local corrections to the Coleman- 
Weinberg potential V e s. 

In figure 2, we plot the Coleman- Weinberg potential V r e ff(</> c ) (resulting from the inte- 
gration of (A.7) for = [i2 = M PL and (0) =V = 3Mp L H 2 / (8tt)) for the same choice 
of parameters. As one can see from figures 1 and 2, the contribution of non-local terms 
to the one-loop effective potential is very small, and hence its effect on the semiclassical 
equation of motion is negligible in a first approximation (see section 4). 

4. Semiclassical field equation and non-local effects 

The analysis of the previous section has shown that the effect of the Bogoliubov coefficients 
on the Coleman- Weinberg potential does not modify the pattern of minima found in static 
de Sitter space [8]. It is therefore legitimate to study the equation (2.12) when, in a 
first approximation, its second line is neglected, and then to use the resulting solution to 
evaluate the non-local correction. Here, the value of V e {{(0) is set equal to Vq as above, 
and it dominates the energy, as it occurs in a de Sitter universe, for 4> c < Mpl- 

On studying the limiting form of equation (2.12) when the second line is neglected, it 
can be easily seen that the self-interaction term (see (A.7)) plays a key role in obtaining a 
sensible physical model. In other words, for vanishing A, equation (2.12) admits a run-away 
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exact solution for C in the form (denoting by Di and D 2 two integration constants) 

C = D 1 exp (-fcHt) + D 2 exp (-(3 2 Ht) (4.1) 

where (3 1 =(3 + V9 + 48Be 2 ) /2 and /3 2 = ^3 - V9 + 48Be 2 ) /2, and we have re-expressed 

r in terms of t by means of (2.10). This behaviour results from the particular form taken 
by the potential in (2.12), i.e. —6BH 2 e 2 (f) 2 , which is unbounded from below. 

We thus study the equation given by the first line of (2.12) with non- vanishing A, i.e. 

1 d 2 6 c 2 da dd> c 9 , . „ , dV e rt , A . 

H — t — ^ — e < A^A^ > C + -^-p- = (4.2) 



a 2 <ir 2 a 3 dr dr c 90, 

subject to the initial conditions C (O) = M PL and d4> c /dr{r = 0) = 0. In figures 3 and 
4, 0c (r) and d4> c /dr are plotted. As one might have expected, c (r) decreases until it 
reaches the zero value, by virtue of the nature of the potential (see figure 2), whilst the 
kinetic energy increases. 

We have then used the solution of equation (4.2) to evaluate the non-local effects 
in our model. For this purpose, starting from the definitions of energy density p^ c and 
pressure p<j> c of the background scalar field in a de Sitter universe, i.e. 



P<f> c = \fc + V eS - 6Be 2 H 2 <p 2 c (4.3) 



P4>c = ^0 2 - Ves + 6Be 2 H 2 tf (4.4) 
one gets from (2.12) the equation 

^ + 3Ha( P<j)c +p+ e ) = -9(0 c (t),t)^. (4.5) 
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With our notation, © corresponds to the whole second line of (2.12), 




(4.6) 



If the right hand side of (4.5) can be viewed as a dissipative term, which implies that it 
always takes negative values and depends quadratically on d4> c /dr, the Bianchi identity 
leads to the following equation for the energy density pr of radiation: 



In figure 5 we plot the right hand side of (4.5) corresponding to the solution of equation 
(4.2). Note that its sign turns out to be positive for the largest part of the r-interval, and 
hence cannot actually lead to dissipative effects. A naturally occurring question is how to 
interpret this lack of dissipation. Indeed, in the adiabatic approximation studied in [5], 
the coupling of the scalar field to a fermionic field by a Yukawa term produces a vacuum 
energy loss rate proportional to (Ha(p c + d(j) c /dT)d4> c /dT. In this last expression the second 
term is clearly a dissipative effect, since it does not depend on the sign of the velocity and 
in the formula reported in [5] it occurs with the correct sign to represent an energy loss. 
By contrast, the first term does depend on the sign of d4> c /dr and in de Sitter, where H 
is constant, it only represents a further quantum correction to the energy and pressure of 
the 4> c fluid. 

These considerations seem to suggest that also in our case, where one deals with decays 
of the classical field configuration into its quanta, the non-local term in the semiclassical 
field equation may be essentially a linear combination of 4> c and d4> c /dr. If this property 



dp R 



+ 4Ha p R = e(Mr),T) 



(4.7) 



dr 



dr ' 
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holds, one can expect that the evaluation of (Hacf) c + d(j) c / dr)d(j) c / dr may indicate when 
dissipative effects are likely to occur, depending on whether the first or the second term of 
this linear combination is dominant. In other words, the sign of the right hand side of (4.5) 
shown in figure 5 can be understood by pointing out that for the solution of (4.2), shown in 
figures 3 and 4, the ratio (dcj) c / dr) / Hacj) c is smaller than 1 whilst d(f> c /dr is negative, when 
r G [0,9.2 Mpl]. Note that, in the neighbourhood of r = 9.2 Mp^, —6 d<p c /dT vanishes, 
and this corresponds to the value of r such that the linear combination (Hacj) c + d(f> c /dr) 
vanishes. 

As a further check of the conjecture about the functional dependence of the non-local 
term of equation (2.12) on 4> c and d(f) c /dT, one should also analyze other situations where, 
unlike before, the ratio (dcpc/dr) / Ha<fi c is larger than 1. This can be done for example by 
taking as a trial function for 4> c (r) the linear combination 4> c (t) = <p c o + ar, where a and 
4> c o are two arbitrary parameters. The resulting analysis, performed by varying a and <p c o 
so as to reproduce the conditions (dcpc/dr) / Ha(p c > 1 or (dcpc/dr) / Ha(p c < 1, seems to 
confirm that, in de Sitter, 0((/> c (r), r) is indeed dominated by a term proportional to the 
combination (Hacj) c + d^ c jdr). 

5. Concluding remarks 

This paper has studied scalar electrodynamics in a spatially flat FRW universe, by in- 
cluding a self-interaction term for a conformally invariant scalar field (cf [5,6,10]). On 
imposing the Lorentz gauge and the supplementary condition (A. 4), the one- loop equation 
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for the expectation value of the complex scalar field in the conformal vacuum shows a 
linear superposition of the self-interaction term of [5] and of the electromagnetic term [6] . 
The results of our investigation are thus as follows. 

First, the numerical solution for the Bogoliubov coefficients in (2.13) and (2.14), sub- 
ject to the initial conditions (2.15), has been obtained and used to find the correction to 
the Coleman- Weinberg effective potential by integrating with respect to <j) c the second line 
of (2.12). The pattern of minima in the effective potential is not modified by the non-local 
term in (2.12). For the particular values of parameters considered in our investigation, the 
non-local corrections turn out to be several orders of magnitude smaller. 

Second, the limiting form of equation (2.12), i.e. equation (4.2), has been studied, 
and its numerical solution has been used to evaluate non-local effects in our cosmological 
model. Such a solution corresponds to a slow-roll dynamics. Interestingly, the approximate 
calculation of the function 6 defined in (4.6) shows that 6 does not lead necessarily to 
dissipative effects in the early universe. Nevertheless, in a de Sitter model, the right hand 
side of equation (4.5) is very well approximated by the same combination of 4> c and d(/) c /dr 
which results from the adiabatic case studied in [5] , where the coupling of a scalar field to 
a fermionic field was instead considered. More precisely, we have found that, for various 
forms of </> c (r), the right hand side of (4.5) is essentially given by 



where A is a positive function of r. The first term in round brackets leads to a further 
quantum correction to the energy density of the scalar field, at least when the function A 
is slowly varying. The second term is purely dissipative (see the end of section 2). 
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One of the main motivations for studying non-local corrections to the Coleman- 
Weinberg potential was their possible relevance for the reheating of the early universe 
[7]. Our analysis confirms that they cannot be simply re-expressed by a term of the kind 
Td(f) c /dT, as first found in [7], where a different approach to the effective action is used 
with respect to [5] . A non-trivial open problem is the numerical evaluation of the function 
A in our de Sitter model. Moreover, it appears necessary to understand how the form of 
A depends on the specific choice of the background field 4> c (t). Last, but not least, the 
whole analysis should be repeated for FRW models which are not in a de Sitter phase. 
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Appendix 

Following [6], we denote by D^{x, x') the photon retarded Green's function, which satisfies 
the equation 

(<$£□ -V»V p -R%) x D%(x,x')=g^5(x,x')/^feU ) . (Al) 
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Thus, by writing A? n (x) for the photon in-field, the solution of (2.4) to order e, in the 
Lorentz gauge 

(A.2) 



= 



is given by [6] 



A^^A^-ie J d 4 x' V-det g(x') D^(x, x') (Vv^ - </>V^ f ) . (A3) 

The residual gauge freedom of the problem is dealt with by imposing the additional con- 
dition [6] 



X»A„ = 



(A4) 



where X is the same timelike vector field appearing in (2.9). Equations (A.2)-(A.4) imply 
that, on inserting (A. 3) into the right-hand side of (2.5), the only non-trivial contribution 
still vanishes after integration by parts, since 



V ^c = V M D% = 



(A5) 



In the equation (2.12) for the expectation value <j) c of the scalar field 0, the derivative 
with respect to cj) c of the one-loop effective potential is given, in the case of non- vanishing 
renormalized mass, by (cf (3.15) of [5]) 

dV eS o , . 1 . 1 . /. 1\_ , 5 



d(p c 



+ 38i^ A ( m2+ ^4) fl + T2^)^ l0g 

+ W A ( m2+ ( ? "5) B + 3^)^ los 



m 2 + «-(l/6))i? + (l/12)A# 



m 2 + (£-(l/6))#+(l/4)A<^ 



(A6) 
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and for m = by (cf (3.16) of [5]) 



dV eS 



d(p c 



4608tt 2 



1 




(g-(l/6))fl+(l/12)A0g _ 11 
(A/12)/*? 3 



) 



+ 5T^ A2 ^('° g 
+ 385? A ( ?_ 5) 

+ 12^ A («-l)- 





(-4-7) 



In agreement with the notation of section 2, m, £ and A are the renormalized values of our 
parameters, and Hi and \ii are completely arbitrary renormalization points. Our equations 
(A. 6) and (A. 7) are obtained by imposing the renormalization conditions (3.12a)-(3.12c) 
of [5] and bearing in mind that the numerical coefficients in our equation (2.8) differ from 
the ones occurring in equation (2.19) of [5], since we study a complex scalar field. 
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Figure captions: 

Figure 1. In Planck units (used also in figures 2,3,4 and 5), the correction to the Coleman- 
Weinberg potential resulting from the Bogoliubov coefficients in (2.12)-(2.14) is plotted 
versus r and 4> c q. 

Figure 2. The Coleman- Weinberg one- loop effective potential is plotted versus 4> c q, in 
the case of a de Sitter background. 

Figure 3. The solution of the semiclassical field equation (4.2) for the background field 
configuration is plotted versus r. 

Figure 4. For the numerical solution plotted in figure 3, its derivative with respect to r 
is here shown. 

Figure 5. The right hand side of equation (4.5) is plotted versus r. 
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